The kagome antiferromagnet: a chiral topological spin liquid? 
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Inspired by the recent discovery of a new instability towards a chiral phase of the classical Heisen- 
berg model on the kagome lattice, we propose a specific chiral spin liquid that reconciles different, 
well-established results concerning both the classical and quantum models. This proposal is analyzed 
in an extended mean-field Schwinger boson framework encompassing time reversal symmetry break- 
ing phases which allows both a classical and a quantum phase description. At low temperatures, 
we find quantum fluctuations favor this chiral phase, which is stable against small perturbations of 
second and third neighbor interactions. For spin-I/2 this phase may be, beyond mean-field, a chiral 
gapped spin liquid. Such a phase is consistent with Density Matrix Renormalization Group results 
of Yan et al. (Science 322, 1173 (2011)). Mysterious features of the low lying excitations of exact 
diagonalization spectra also find an explanation in this framework. Moreover, thermal fluctuations 
compete with quantum ones and induce a transition from this flux phase to a planar zero flux phase 
at a non zero value of the renormalized temperature (T/S^), reconciling these results with those 
obtained for the classical system. 



PACS numbers: 75.50.Ee,75.10.Hk,75.40.Cx 

With the extensive degeneracy of its classical ground- 
state, the antiferromagnetic Heisenberg model on the 
kagome lattice was recognized early on as the paradigm 
of a quantum spin liquid phase. ' The recent experimen- 
tal discovery of Hcrbcrsmithite, a kagome compound 
fluctuating down to temperatures thousands of times 
lower than the coupling constant, strengthens this spec- 
ulation.' Despite numerous efforts, the nature of the 
ground state (GS) of the spin-1/2 Heisenberg antiferro- 
magnetic kagome model (AFKM) remains controversial. 
Exact quantum approaches point to the absence of long 
range order. ' Although exact diagonalizations (ED) on 
small samples (up to N=36) leave open the question of 
the criticality''' ' , Density Matrix Renormalization Group 
(DMRG) calculations^ support the idea of a true gapped 
spin liquid. 

Recently, a new instability of the degenerate classical 
model towards a chiral phase has been discovered. In 
this paper, we show at a mean-field level that the hypoth- 
esis of a chiral spin liquid holds and is consistent with nu- 
merous robust results accumulated during the last twenty 
years, both for spin-1/2 and in the classical limit. 

The properties of chiral spin states, with simultane- 
ously and spontaneously broken space reflection (P) and 
time reversal (T) symmetries, were largely debated at the 
end of the eighties in the wake of the Quantum Hall effect. 
A revival of these topics has occured thanks to graphene 
and flat band insulators. Wen et al. ' ' deflned the chiral 
phases through the fluxes of the underlying gauge fields 
and Kalmeyer and Laughlin ' ' proposed to describe spin 
liquids by Laughlin wave functions. Yang et al. sug- 
gested that the Heisenberg model on the kagome lattice 
might be in a chiral spin liquid state. We reexamine 
this suggestion inspired by the knowledge of the classical 
non-planar spin order, described in Fig. 1, and propose 
a specific chiral spin liquid as the GS of the spin-1/2 




FIG. 1: (Color online) Description of the cubocl order. Left; 
on the kagome lattice each color corresponds to a different 
magnetic sublattice. The thick line indicates the 12 site unit- 
cell. Right: arrows are the spin orientations with the same 
color coding as in the left flgure. The black lines connecting 
the end of the vectors form a cuboctaedron. On each triangle, 
the spins are coplanar at 120 degrees; for opposite sites on 
each hexagon, spins are anti-parallel. The triple products 
(determinant) of three spins of the hexagons, either flrst or 
second neighbors, are non zero and measure the chirality of 
the phase. They change sign in a mirror symmetry or in a 
spin flip. 

AFKM. 

The kagome lattice can be viewed as a lattice of corner 
sharing triangles. The classical ground-state on a single 
triangle is planar with three spins at 120 degrees. Fixing 
the spin plane on a triangle does not fix the planes on 
adjacent ones whence the extensive ground-state degen- 
eracy. This degeneracy can be lifted by couplings beyond 
nearest neighbors. We consider the following Ji-J2-Jzh 
Hamiltonian with Ji = 1: 

H=Y. jj,^ ■ S, (1) 
where the coupling constants and exact classical phase 
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diagram are given in Fig. 2. AFKM refers to this model 
when J2 = Ja/i = 0. For J^h = 0, numerous studies of 
quantum states have been motivated by the two planar 
classical states denoted q = (J2 > 0) and \/3 x -\/3 
(J2 < 0)^'^''^"''^'"^''. An infinitesimal antiferromagnetic 
third-neighbor interaction across hexagons ( J3/,. > 0) lifts 
the degeneracy of the classical AFKM to a 12-sublattice 
magnetic state where the spins point toward the corners 
of a cuboctaedron (Fig. 1), whence the name cuhocl. 
This order was first introduced by Janson et al. "'" who 
claim that this J^h > interaction should be of experi- 
mental relevance. This order is chiral: it breaks mirror 
symmetry. 

Monitoring the evolution under the effect of quantum 
fluctuations from the classical limit to the disordered 
spin-1/2 system remains a challenge. The Schwinger bo- 
son mean-field theory (SBMFT) is an approximate but 
versatile method to study, in an unified framework, both 
long range ordered (LRO) and gapped spin liquid phases, 
from classical to the quantum limit. In a first enlighten- 
ing work, Sachdev ' ' showed that amongst the planar 
states, the y/3 x VS is more stable than the q = at 
the AFKM point. In order to study the present model 
(Eq. 1) around the AFKM point, we extend Sachdev's 
work in several directions as will be seen below. 

The Schwinger boson operator bl^ {a =t or l) creates 
a spin 1/2 on lattice site i. A physical spin S at site i is 
represented by 2S bosons. After a mean field decoupling, 
the Hamiltonian reads: 

Hmf = Jc,{BijB\^ - AjA\j) + h.c. - ^ XiTi, + eo, 

. . (2) 

where the bond operators are defined by 2Aij = bi^bji^ — 

bij)j^, 2Bij = bl^bj^ + blj)ji- Atj and Bij are the 
associated complex mean-field parameters to be deter- 
mined by the self-consistency equations. The {Aij,Bij} 
set is called an ansatz. are Lagrange multipliers 
to constrain the mean boson number: (n^) ~ 2S and 
£0 = Jai\Aij\'^ - \Bij\'^) + 2S J2t Aj, where 5 is a 

continuous real positive mean field parameter. 

Most SBMFT studies use only one of the two types of 
parameters {A or B). Recently, taking both fields, Mezio 
et al. found a much better description of the excitation 
spectrum of frustrated systems. More specifically, the 
A fields describe the singlet amplitudes whereas the B 
fields allow the description of boson hopping amplitudes 
that are a fundamental ingredient to describe the mixing 
of spin singlets and triplets on each bond, a mechanism 
that is central in quantum frustrated magnets. In addi- 
tion, because both ferromagnetic and antiferromagnetic 
interactions are treated on an equal footing, the phase 
diagram can be explored continuously around the origin 
regardless of the sign of coupling parameters. 

Solving the full problem with two complex parameters 
per link and one real Lagrange multiplier per site is nu- 
merically too demanding for large lattices. Looking for 
spin liquids or regular LRO we assume = A and 
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FIG. 2; (Color online) Phase diagram of the model with 
up to third neighbors interactions at T = and Ji = 1. (a) 
exact classical phase diagram (the method used to obtain this 
phase diagram and all orders are described in ). The point 
AFKM: (J2, Jsh) = (0,0) is a tricritical point. (6) SBMFT 
phase diagram for S = 0.5. The tricritical point stands at 
{J2,J'ih) = (0.0049, -0.021) and the AFKM point (red circle) 
is now inside the cubocl phase. 

A and B are invariant under lattice symmetries up to a 
local gauge transformation. Using projective symmetry 
groups (PSG) , Wang and Vishwanath ' ^ obtained four 
ansatzes where physical observables are invariant under 
lattice symmetries. They are defined by the fiuxes of the 
A operators on specific loops ((/)o,(/>^) = (0,0), (7r,0) , 
(0, tt) and (tt, tt). y/S x -\/3 and q = are associated with 
the first two choices respectively. 

The chiral cubocl ansatz, however, cannot be obtained 
within this first PSG approach, because in chiral states 
the lattice symmetries are respected only up to a time 
reversal symmetry. In all previous studies, the mean- 
field parameters were chosen real, the fiuxes equal to 
or TT, thus excluding chiral ansatzes. The extension of the 
PSG to include both the symmetric and the chiral spin 
liquids will be described in a longer paper-^^ . In short, the 
new ansatzes are defined by complex fields with specific 
constraints on the moduli and arguments. Thanks to the 
PSG analysis the number of parameters at the AFKM 
point is limited to 2 moduli of bond fields for q = or 
x -\/3 plus a phase Oai^ ^ot cubocl. Other bond fields 
of the 6-spin unit cell are fixed by algebraic constraints 
(see Fig. 3 and Table I). Non zero fluxes ' (modulo tt) 
induced by ^Ai^ are an indirect mean-field measure of the 
chirality of cubocl, (they are or vr for q = or \/3 x \/3 
respectively) . 

The numerical solution of the mean-field equation is 
found from a descent method minimizing the sum S of 
the squares of the (free) energy derivatives with respect 
to the field parameters, each single energy evaluation be- 
ing maximized with respect to A (we stop the descent 
when E < 10~*). We keep the solutions with hessians 



FIG. 3: (Color online) Unit-cell (in thin green lines) of the 
cubocl ansatz (left) and q = and \/3 x \/3 ansatzes (right) 
at the AFKM point. Aia, An,, Bia and Bib are complex link 
parameters with constraints on their moduli; |j4iaj = 
\Bia\ = \Bii,\ and on their arguments ~ 0, Obi^ ~ = 
TT. The last constraint depend on the ansatz: Oa^^ = for 
q = , TT for \/3 x \/3 and is not fixed for cubocl. 
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TABLE I: Values of the self-consistent SBMFT parameters 
for the three competing ansatzes near the AFKM point. "I^The 
exact value is tt — arctan y/S. 



of the correct sign, positive for the A fields and negative 
for the S if J > 0, and the opposite if J < 0. These 
requirements imply that our solutions are stable against 
gaussian fluctuations. 

The resulting phase diagram at T = and S = 0.5 
is given in Fig. 2-b: at the AFKM point, the cubocl 
ansatz is more stable than any other regular ansatz, with 
an energy per site of —0.4717. The numerical values 
of the parameters at this point are given in Tab. I. The 
parameter range of stability of the cubocl phase increases 
when the spin decreases from S = 0.5 to 0.366. For J^^ — 
and 5 = 0.5 it is Ja G [-0.005; 0.025] and for S = 0.366 
it is enlarged to J2 £ [—0.008; 0.045]. This increase is 
another proof of the role of quantum fluctuations in the 
stabilization of the cubocl phase. 

The dimensionless free energy difference (AF/iS^) be- 
tween the \/3 X a/3 phase and the cubocl phase is given 
as a function of S in Fig. 4-(a). At T = 0, it is of the 
order of 10""^ in favor of cubocl. In the classical limit, 
5 — cx), the two phases are degenerate as they should 
be. The comparison with the q = state is not shown as 
it always has a much higher energy at the AFKM point. 

Decreasing S leads to a second order phase transition 
from a gapless LRO cubocl phase to a fully gapped chiral 



FIG. 4: (a) Difference of SBMFT free energy between the 
cubocl and \/3 x \/3 ansatzes as a function of 5 at T = 0. 
Inset: zoom of the domain around critical spins 5c. (b) Same 
quantity as a function of T for different 5. Inset: value of 
Tc/S'^ versus 5 (see the text). 



spin liquid at a critical value Sc = 0.4. One should not 
hastily conclude that the true spin- 1/2 system has Neel 
long range order. In this mean-field approach, the on-site 
number of bosons fluctuates: it is only fixed on average, 
5 is a parameter and (S^) = 3S{S + l)/2 . For spin- 
1/2, the good quantum number is (S^) =3/4. To recover 
this good quantum number we should use the parameter 
5 = (V3 - l)/2 ~ 0.366. With S = 0.5 the phase is 
gapless and finite size-scaling shows a very small stiffness, 
while with S 0.366 the system is a gapped spin liquid 
compatible with the results of Yan et al. . 

Now, we turn to the effect of thermal fiuctuations at 
the AFKM point. Fig. 4-b shows how they destabilize 
the cubocl phase in favor of the y/2> x one. The renor- 
malized transition temperature Tc/S"^ from the cubocl to 
a/3 X a/3 phase decreases to zero with increasing S (Inset 
Fig. 4-b). This is consistent with classical numerical sim- 
ulations showing a selection of a planar state by thermal 
fiuctuations. 

This chiral order hypothesis also explains why the ED 
spectra for sizes up to A'' = 36 have a large number of 
singlets below the triplet gap. Let us consider a cell of 
12 spin-1/2 describing some short range order (SRO). As- 
suming three spin directions (as for y/i x a/3 and q = 
SRO), one finds a single S=Q ground-state derived from 
the coupling of three spin-2 dressed by quantum fluc- 
tuations. With a 12-sublattice cubocl SRO, one finds 
132 singlets built starting from the angular addition of 
12 spin-1/2, many of which are low energy states. This 
crude picture makes it possible to understand why there 
are so many (63) singlet states below the triplet gap of 
the A'^ = 36 sample. Moreover this property may explain, 
through resonances, the stabilization of cubocl relative 
to V3 X V3 SRO. 

On the other hand with the hypothesis of a chiral spin 
liquid we expect an 8-fold GS degeneracy on a 2-torus (a 
factor 2 for the chirality times a factor 4 for the topo- 
logical degeneracy) at the thermodynamic limit. Thus, 
the large number of low lying singlets seen in ED spec- 
tra should be restricted to small size samples. Such an 
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evolution has already been observed in the J1-J2 model 
on the triangular lattice in the parameter range where 
the classical ground-state has a 4-sublattice unit cell. 
Note that DMRG results do not exhibit a large number 
of singlet states below the triplet gap.'' 

Spin- 1/2 ED results on the = 36 sample are also 
compatible with cubocl short range order: i) the first S = 
1-state is at the softest fc-vector of the cubocl short range 
order as can be seen in Fig.4 of ref. '; ii) the dynamical 
and static structure factors"'^ have relatively larger values 
at the wave vector of the cubocl order than near the quasi 
soft points of the ^/3 x ^/3 and q = orders. Thus, 
correlation functions in large scale DMRG computations 
and/or characterization of low energy excitations by ED 
for 48 sites samples would be an essential complement to 
further support or discard the present proposal. 

Moreover, the first spin-1/2 states of small samples 
(with an odd-integer number of sites) have non-zero 
Chern numbers. ' This quantum number, first introduced 
in such a context by Haldane and Arovas"'', is a topo- 
logical index (and thus a robust property) characteriz- 
ing the chiral character of a wave function. This prop- 



erty, which has never been explained in other approaches, 
could be understood for chiral spinons.' ' In the same 
spirit, the classical chirality defined by the determinant 
of three spins (non zero on hexagons for the classical 
cubocl ) is generalized for quantum spins as the imagi- 
nary part of the cyclic permutation operator of spins on 
closed contours. This quantity (Wilson loop operator) 
computed in the ED-GS for different contours obeys the 
law expected for a chiral liquid"' ' . 

We have shown, that within mean field Schwinger bo- 
son approximation, the kagome antiferromagnet has a 
chiral ground state. Spin-1/2 exact results both from ED 
and DMRG give some support to this hypothesis. In such 
a system, we expect a low temperature chiral symmetry 
breaking phase with topologically protected edge states, 
as in the quantum Hall systems. The extent of this low 
temperature phase, the nature of the phase transition, 
and the role of defects,- questions already addressed in 
classical systems remain open questions. 

We thank Gregoire Misguich and Frederic Mila for 
many stimulating discussions. We acknowledge Julian 
Talbot for a critical reading of the manuscript. 
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